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Abstract 

In experiments which take a hquid crystal rapidly from the isotropic to the ne- 
matic phase, a dense tangle of defects is formed. In nematics, there are in principle 
both line and point defects ("hedgehogs"), but no point defects are observed until 
the defect network has coarsened appreciably. In this letter the expected density 
of point defects is shown to be extremely low, approximately 10~^ per initially cor- 
related domain, as result of the topology (specifically, the homology) of the order 
parameter space. 
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An outstanding puzzle in the formation of defects after rapid quenches is the ab- 
sence of point dischnations ("hedgehogs" or "monopoles" ) in nematic hquid crystals. 
One might naively expect of order one defect per initially correlated domain. However, 
experiments with rapid pressure quenches do not find any monopoles at all to begin 
with, although the expected dense network of line dischnations ("strings") is present, and 
rapidly reaches a scaling regime in which the string length density decreases as t~^. This 
can be understood to be a result of the t^^^ growth law of the network scale length L{t). 
If there is on average about 1 segment of string of length L per volume L^, the scaling of 
the length density is explained. The monopoles do not make their appearance until 
relatively late in the scaling regime, forming from collapsing loops of string. 

In this letter a solution to this puzzle is presented. The answer lies in the topology of 
the order parameter space, which is the projective plane RP^. In order for there to be a 
monopole inside some sphere in the liquid crystal, the order parameter field has to cover 
its entire space twice. It turns out that this is very hard to arrange out of the random 
initial conditions produced by a rapid quench. There is an underlying mathematical 
formulation of the solution, in terms of the homology of the order parameter space, 
which is outlined briefly at the end of this work. 

The order parameter of a liquid crystal is a traceless symmetric rank 2 tensor Qij{x). 
The normalized eigenvector with the largest eigenvalue is known as the director field 
ni{x), for it defines the average local orientation of the liquid crystal molecules. In a 
nematic the other two eigenvalues are equal, and we may write 

Qij = A^niUj -^6ij). (1) 

The free energy of the liquid crystal is 

F[Q] = J dPx {LidkQijdkQij + L2djQijdkQik + L^dkQijdjQik + V{Q)) (2) 

where V{Q) is the bulk free energy. Near the phase transition we are justified in expand- 
ing to quartic order, and 

ViQ) = \a tr (g^) + 1/3 tr (Q^) + I7 ^AQ'f + ■■■■ (3) 

The condition that the system be in the nematic phase, i.e. that A ^ minimize the free 
energy, is just a'y/ (3"^ < 1/9. In this phase the symmetry group of the bulk free energy 
density, which is the group of spatial rotations SO (3), is reduced to the cylinder group 
DoQ, or 0(2). The manifold M of possible equilibrium states is defined by the condition 
6F/6Q = 0, subject to the constraints of tracelessness and symmetry. This is isomorphic 
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to the coset space SO(3)/0(2), or the real projective plane RP^, which can be thought 
of as a 2-sphere with antipodal points identified. 

After a rapid quench, the order parameter is uncorrelated beyond a certain distance 
.^0, which is determined by the relative magnitudes of the quench time and the relaxation 
time of the system. The isotropic-nematic transition is weakly first order, which means 
that the correlation length grows by a large factor as the phase transition is approached 
0], although the transition itself appears to proceed by bubble nucleation and growth p. 
The substance of the present argument should not depend on the order of the transition: 
it is only the existence of uncorrelated domains that is important. Defects form at the 
interstices of domains where the order parameter is in some sense maximally misaligned, 
in what is known in cosmology as the Kibble mechanism One can estimate the density 
of defects by applying the so-called "geodesic rule" @, |, This assumes that if we pick 
two points Xi and X2 in adjacent domains with order parameters and Q^, the most 
likely interpolation Q^'^{x) on a line between the points is the shortest path in M, for this 
locally minimizes the bulk free energy. Since M comes equipped with a metric by virtue 
of its embedding in the Euclidean field space, this path is by definition a geodesic. Now 
consider three adjacent domains, and pick three points {a^i, X2, x^}. The geodesic rule can 
be applied separately to each pair of domains, and then to all three: the interpolation 
Q^'^^{x) to the interior of the triangle {xi,X2.,x^} is a geodesic surface in M. There 
may be an obstruction to the procedure: if the loop {Q^^ ■, Q^^, Q^^} is in the non-trivial 
homotopy class of 7ri(M), a line defect must pass through the triangle {a^i, X2, x^}, at the 
junction of the three correlated domains. A similar argument involving four domains is 
applied to the formation of point defects 0] , which are obstructions to the construction of 
an interpolating geodesic 3-simplex Q^^^^. We shall see, however, that four uncorrelated 
domains are not enough for a point disclination in a nematic liquid crystal. 

Calculating the probability of finding a defect associated with non-trivial 7r„(M) at 
the interstices of n + 1 domains is a problem in geometric probability on the manifold 
of equilibrium states M. This problem has been solved only for M and for 1- 

dimensional defects in RP^ and S'^ /Z2 The solution is rather neat for the spheres. 
Consider first = 1, where the order parameter is a 2 component field (pa with 0^ 
constant. The problem consists essentially of placing 3 points 0^, 0^, 0^ at random on 
the circle of constant I^a^a; ^^"^ asking the probability for (fy^ to lie between —0^ and 
— 0^ (taking the shortest route). In that case, and in only that case, will the geodesic rule 
supply a loop which wraps around M. Now, ±0^ and ±0^ divide the circle into 4. Given 
that 0^ and cfP' are isotropically distributed, one can convince oneself that the average 
length of the line segment between —0^ and —0^ is 1/4. This is then the probability of 
finding a line defect at the junction of three adjacent domains, and the number of defects 



3 



per unit area is therefore 1/4^q. This generahzes for arbitrary n to 1/2"+^. For strings 
in RP^ the calculation is more involved, but it emerges that the probability is l/vr. 

The problem with trying to extend these calculations to point defects in RP^ is that 
four neighbouring uncorrelated domains can never generate such a defect. To construct 
a hedgehog configuration of the order parameter we must cover M twice, because the 
director field has an x —x symmetry. One cannot unambiguously do this with four 
domains, for the geodesic rule produces a mapping from the tetrahedron {a^i, 0^3, 3^4} 
which is either trivial or contains a string passing through two of the faces. The point is 
that in order to cover M twice, each face of the tetrahedron must cover on average half 
of it, which means that there will always be faces trying to cover more than half. This 
cannot happen with the geodesic rule. Thus we need more domains, which inevitably 
lowers the probability of finding a defect. 

The minimal triangulation of RP^ has in fact 6 vertices (see Figure 1). One can think 
of this as a triangulation of S"^ by an icosahedron, with antipodal points then identified. 
Thus in order to cover RP'^ twice we need a roughly spherical arrangement of a minimum 
of 12 uncorrelated adjacent domains. A great deal of calculation can now be saved by 
an approximation which uses a fixed triangulation of RP^ directly. For example |T0| , pA| , 
if we approximate by three equidistant points labelled 0,1,2, and assign a string to 
a spatial triangle {a^i, a;2, a^s} when all three values of are different, the probability of 
having a string passing through the triangle is just the number of different arrangements 
of 0,1, and 2 divided by the total number of possible assignments 3^. Thus the probability 
in this discrete approximation of 1-dimensional defect passing through the triangle is 

P[[S^) = 3!/3^ = 2/9, (4) 

where the prime is used to denote the approximation to the true geometric probabilities 
Pn{M). For general n we have 

P^(5") = (n + 2)!/(n + 2ri (5) 

This approximation gets worse for large n. Using Stirling's approximation, one sees that 

For line defects in RP^ the calculation proceeds as follows. The first two values of the 
order parameter and can be any two different vertices of the triangulation. The 
last point must be one of the two which are connected to both of the first two. Thus 

P[{RP^) = 6 ■ 5 ■ 2/6^ = 10/36 (6) 

which is close to the true value Pi{RP^) = l/vr. For point defects, we must calculate the 
number of different ways of assigning values of Q to the 12 domains. Picking any two 
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adjacent domains, the first values can once again be any two vertices of tlie triangulation. 
In a tliird domain, adjacent to botli tlie first two, one must correspondingly pick one of 
the vertices connected to both those already selected. Thus 

P^{RP^) = 60/6^^ ~ 2.76 X 10"^ (7) 

A quick way of calculating this number is to note that the assignment of vertices to 
domains is just a map from one icosahedron to another with opposite points identified. 
Therefore P2{RP^) is just the order of the icosahedral group, which is 120, divided by 2. 

The configuration of domains occupies a volume of approximately ^q, and so the 
density of point defects Np is roughly 

Np ^ 10-%-l (8) 

This is a very small number, as promised. If the discrete calculation is here as good 
an approximation as for the spheres, then it explains why the point defects of a nematic 
liquid crystal are not found after a rapid quench: they require a very special arrangement 
of the order parameter over many uncorrelated domains W2 . 



The icosahedral arrangement of domains can be extended into the body of the material 
by the addition of a further domain in the centre. One then realises that the "point" 
disclination is actually a small loop of size ~ encircling the central domain. The value 
of the order parameter here merely controls the loop's orientation. Thus there is a sense 
in which there are no point disclinations at all. What we have calculated is merely the 
density of the smallest possible loops which can form hedgehogs. 

To conclude, I outline the mathematical structure implicit in the geodesic rule. Recall 
that the construction starts with points {xi] in uncorrelated domains, and the corre- 
sponding values of the order parameter {Q*}. One attempts to construct an approxi- 
mation to the field configuration over the whole of by extending the points {xi] to 
a full triangulation, defining the order parameter field Q{x) by the geodesic rule. This 
determines how to "fill in" the set of closed figures (points, lines, and triangles) in order 
to create others (lines, triangles and tetrahedra) of higher dimension. The result is a 
simplicial complex |]l3l in the order parameter space M. However, the procedure fails 
when some sub-complex cannot be filled in, that is, the sub-complex is not the boundary 
of another, higher- dimensional complex in the space M. The order parameter has to 
leave M and a defect appears in the corresponding region of R^. This can happen if and 
only if the space has a non-trivial homology group Hn{M). Thus the Kibble mechanism 
coupled with the geodesic rule produces defects of dimension ci in a space of dimension 
D only if Hn{M), with n = D — d—1, is non-trivial. The second homology class of RP^ 
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is zero: this is the underlying reason for the absence of point defects in nematic hquid 
crystals. 

I am indebted to Neil Turok for many useful discussions. I also thank Martin Za- 
potocky for encouragement to write this work up, and for a careful reading of the first 
draft of the manuscript. I also wish to acknowledge the hospitality of the Isaac Newton 
Institute where this paper was completed. This work is supported by PPARC Advanced 
Research Fellowship B/93/AF/1642. 
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Figure Captions 

Figure 1: The minimal triangulation of RP^, consisting of 6 vertices, 12 edges and 10 
faces. This is essentially the top half of an icosahedron. Opposite points on the boundary 
are identified. 
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Figure 1 



